Genome Biology and Evolution

BackgroundAlthough strict maternal transmission of mitochondria is a general feature of animals and humans for ensuring homogeneity in mitochondrial DNA (mtDNA) across generations, exceptions were reported in the recent past. For example, some extremely rare but spectacular cases of heteroplasmy and paternal transmission in humans have questioned the universal evolutionary principle. Hence, as an alternative, the Mega-NUMT concept was coined to explain this discovery and was thereafter partly proven to exist. This concept expands on the quite common transfer of mtDNA fragments to the nucleus (NUMTs) by considering the existence of multicopy mitochondrial nuclear insertions. Mega-NUMT reports are currently restricted to a few cases in animals, including humans. However, even in humans, their detailed genomic organization, natural prevalence, and potential biological functions remain unclear. Methodology/Principal FindingsHere, we discovered that up to 60 full-sized mitochondrial genomes are integrated into the nuclear genome of the neotropical fruit fly Drosophila paulistorum using long-read sequencing and confirmed their presence by in situ hybridization. The copies are organized in one cluster on chromosome 3, which we, due to its similarity with the Mega-NUMT concept, designated the "Dpau Mega-NUMT". Contrary to the rarity in humans, this Mega-NUMT is found at high prevalence (40%) in both long-term laboratory lines and natural D. paulistorum populations of different semispecies. Additionally, the mitochondrial copies in the Mega-NUMT cluster are phylogenetically separated from the current mitotypes of D. paulistorum. Together, these observations suggest long-term maintenance of the Mega-NUMT in nature. Hence, we propose that the Dpau Mega-NUMT may have been transferred to the nuclear genome before D. paulistorum semispecies radiation and maintained at relatively high prevalence in nature by balancing selection due to yet undetermined functions. Conclusions/SignificanceTo our knowledge, this is the first verified existence and detailed dissection of a Mega-NUMT outside cats and humans. We show that Mega-NUMTs can be persistent in nature, even at high prevalence, potentially due to balancing selection. Our findings strengthen the importance of high-quality long-read sequencing technologies for deciphering complex repeat-rich genomic regions to deepen our understanding of the dynamics of genome evolution within genomic "dark matter".

Successful establishment of long-term, obligate endosymbiotic relationships requires integration of hosts and endosymbionts across multiple levels. For example, highly integrated, host-beneficial endosymbionts typically have extremely reduced genomes and metabolisms. However, we do not yet fully understand the specific mechanisms that drive this integration or if there is a specific order in which these changes must occur. To investigate the early stages of endosymbiont genome reduction, we greatly expanded available whole genome data for the nitrogen-fixing endosymbionts (spheroid bodies, SBs) of diatoms in the family Rhopalodiaceae. We used these data to reconstruct SB evolutionary history and to characterize SB core metabolic capacity. We found two key genes missing from all SB genomes, mltA and dnaA, which could provide points of host control over SB cell division. Although most of the SB core genome is experiencing moderately strong purifying selection, we identified 54 genes under positive selection. Eighteen of these are peripheral proteins or involved in cell wall and cell membrane metabolism and could be involved in direct interactions with the host. Unexpectedly, we also found three nif genes under positive selection that are core to the central nitrogen-fixing enzyme. Overall, our results provide early insights into how SBs and their hosts interact, showing that SBs are still in the early stages of endosymbiont genome reduction, but they differ in key ways from current models, including the early loss of all mobile elements.

Population genetic theory predicts that natural selection will be more efficient in large than small populations because genetic drift reduces the efficiency of selection in small populations. Small populations adapting to new environments can also be expected to evolve higher recombination rates to facilitate adaptation as well as to dissociate and purge harmful mutations. We tested these hypotheses (1) by investigating differences in the strength of association between nucleotide diversity ({pi}) and recombination rate across the genomes of nine-spined sticklebacks (Pungitius pungitius) from four small freshwater (mean Ne {approx} 2 578) and four large marine (mean Ne = 86 742) populations, as well as (2) by comparing recombination rates between small and large populations using population specific linkage maps. We found the predicted positive correlation of{pi} with recombination rate from all but the smallest freshwater populations, suggesting prevalent linked selection even after accounting for variation in GC/CpG content, and gene density. Mean recombination rates did not differ between freshwater and marine populations, except that the smallest Ne freshwater population exhibited significantly elevated recombination rate. GWAS analyses suggested a polygenic basis for recombination rates. These results suggest an important role for linked selection in reducing{pi} in low recombination regions especially in large populations. Moreover, as predicted by theory, at least one of the small freshwater populations appears to have evolved a higher recombination rate than its marine ancestors.

Accurate estimation of population demographic history is central to population genetics yet remains challenging due to the sensitivity of inference methods to the number of individuals and the demographic scenario assumed in inference. The site-frequency spectrum (SFS) of neutral variants, a widely used summary statistic of genetic variation, is particularly sensitive to demographic processes, but studies have shown that qualitative results from demographic inference, i.e., population expansion vs. contraction, can depend strongly on the number of individuals in the dataset. Here, we analyzed two simulated datasets and one empirical dataset characterized by an ancient population bottleneck followed by a recent population expansion. Fitting a two-epoch demographic model across a range of sample sizes, we found that inference shifted from signals of ancient population contraction at small sample sizes to signals of recent population expansion at large sample sizes. Other summary statistics, including Tajimas D and the proportion of singletons, also changed with sample size. We found that these changes of inferred evolutionary signals under a two-epoch model can be explained by the epoch which contributes the highest mean proportion of coalescent branch lengths. Our results highlight that demographic inference depends critically on the number of individuals analyzed and suggest that analyzing datasets at multiple sample sizes can reveal complementary aspects of population history.

Small marsupials in the family Dasyuridae are a key component of Australias arid and semi-arid fauna, whose high species richness is proposed to reflect an opportunity-driven adaptive radiation. Despite growing interest in this group from both ecological and evolutionary perspectives, genomic data for most species is non-existent, or limited to a few marker loci. Here, we generated a chromosome-level reference genome and a de novo mitochondrial genome for the desert-dwelling Wongai ningaui (Ningaui ridei). The nuclear genome assembly is highly contiguous, with a scaffold N50 of 594.484 MB and high BUSCO gene recovery (93.84%). Additionally, we produced a draft assembly for the related, semi-arid slender-tailed dunnart (Sminthopsis murina). We then used these assemblies to explore the demographic histories of these species. We find evidence for contrasting patterns of population growth during the late Pleistocene and early Holocene, corresponding with differences in local climate, potentially consistent with differences in optimal habitat. The new genomic resources and demographic findings presented here provide a foundation for future studies on adaptive specialisation in this group of Australian marsupials. Significance StatementDasyurid marsupials are the primary carnivorous and insectivorous mammals in Australia. This diverse family includes species such as the endangered Tasmanian devil (Sarcophilus harrisii) and quolls (Genus Dasyurus), as well as an emerging laboratory model species, the fat-tailed dunnart (Sminthopsis crassicaudata). Despite the species richness within dasyurids, most species remain under-studied. This is particularly true of arid and semi-arid zone species, who are often small in size, live in remote habitats and are cryptic by nature. By creating genome assemblies for two dasyurid species, this study provides resources to support a variety of phylogenetic, population genetic and evolutionary developmental lines of research. Importantly, the studys finding that arid and semi-arid dasyurids show distinct trajectories of demographic change in response to historical climatic shifts may point to local adaptations with implications for the resilience of these species to ongoing and future climate change.

Global populations are ageing at an unprecedented rate. For many diseases, age is a strong indicator of susceptibility, morbidity, or mortality. Principles of evolutionary biology can be leveraged to understand how pathogens may optimally exploit new populations, and the impact of this on the global burden of infectious-disease-induced mortality. We parameterised an age-specific R0 model with 2017 epidemiological data on Measles, Tuberculosis, Meningitis, and Ebola, and age-specific demographic estimates for 2017 and 2050, for the seven Global Burden of Disease super-regions. We explored the theoretical trade-offs between pathogen virulence & transmission, and virulence & host recovery, parameterising trade-off parameters using Latin Hypercube Sampling. Population ageing between 2017-2050 saw an increase in virulence induced mortality in four settings: 1) Ebola in sub-Saharan Africa, 2) Measles in central/eastern Europe & central Asia region, 3) Measles in North Africa & the Middle East and 4) Tuberculosis in the central/eastern Europe & central Asia region. The decrease in infection duration due to an increase of elderly people drives pathogen virulence down for diseases in the remaining settings. Understanding the mechanisms that shape pathogen dynamics and leveraging this to predict future challenges is key to endemic disease management in a rapidly changing world. Author SummaryKey aspects of disease transmission including susceptibility to infection, the severity of infection, and the probability of dying from that infection, are affected by host age. Global populations are rapidly ageing, so that the mean age of most populations is generally higher than it used to be and is set to continue on this trajectory. This suggests that the dynamics of infectious diseases are also likely to change, although infectious disease dynamics tend to be non-linear as these key parameters interact. We have developed a dynamic modelling framework to explore how changes in population age structure might impact the optimal level of pathogen virulence in a population. We have chosen four infectious diseases as case studies, that differentially impact certain age classes to illustrate these dynamics. We have parameterised this framework with open access data for each of the seven Global Burden of Disease super-regions and show that population ageing can increase virulence for several diseases in differing global regions, whilst increased background rates of mortality can drive virulence down in others.

Life history traits are often correlated, creating trade-offs that may impede the response to natural selection and be responsible for the evolution of senescence. These trade-offs may arise through pleiotropic effects, which can affect the response to selection in ways that resemble intra-locus sexual antagonism. Despite these hypothesized relationships, we lack clear connections between pleiotropy, sexual antagonism, and the evolution of life histories. Empirical tests for inter-sexual differences in life-history traits, including sex-specific aging, can be used to evaluate hypotheses about how pleiotropy and sexual conflict affect evolutionary trade-offs. To those ends, we measured lifespan, development time, and body size in Drosophila pseudoobscura males and females, each of which carried one of six third chromosome inversion genotypes. Temperature affected lifespan and development more than any other factor; higher temperatures increased mortality rate, decreased lifespan, and accelerated development. However, we also observed sex differences in mortality rates and development times that depended on genotype and temperature. Notably, temperature elevated the initial mortality rate across all flies, yet increasing temperatures reduced the rate of aging in some genotype-sex combinations. Similarly, direct effects of genotype on mortality rate and development time depended greatly on sex and temperature, but there was no genotype effect on body size. Despite these context-dependent genotype effects on life history traits, we failed to identify any correlations that would serve as clear evidence for sexual conflict or trade-offs. Our results therefore suggest that either historical conflicts have been resolved or any conflicts that may exist do not result in the correlations predicted by existing models.

Genetic evidence for fluctuating selection has begun to accumulate for different species over the past few decades, especially for the Drosophila genus where studies have reported hundreds of loci undergoing putatively adaptive oscillations across successive seasons. However, most theoretical and simulation studies of fluctuating selection have relied on abstract or weakly parameterized models, making it difficult to assess their relevance for natural populations. In this study, we simulate multilocus seasonally fluctuating selection under a recently developed model and examine its effect on the variance effective population size (Ne) at a genome-wide scale. By recapitulating genomic, demographic, and evolutionary parameters from natural Drosophila populations in our simulations, we were able to reproduce allele frequency oscillations reported in recent studies and show that these lead to [~]50% genome-wide reductions in Ne. We also demonstrate that Ne reductions are well predicted by the maximum frequency amplitude among all adaptively fluctuating loci, and that the frequency amplitudes are largely determined by the number of adaptively fluctuating loci and the strength of their epistatic interactions. Our results demonstrate that fluctuating selection can substantially reduce effective population size and underscore the importance of temporally variable selection in shaping genome-wide patterns of variation beyond classical models. Article SummaryGenetic studies of fluctuating selection in natural populations have grown steadily over the past decade, with reports suggesting that hundreds of loci undergo adaptive oscillations over seasonal timescales in cosmopolitan Drosophila populations. By simulating seasonally fluctuating selection under a recently developed model and ecological scenarios informed by published studies, the authors show that this mode of selection can reduce effective population size by [~]50%, with the magnitude of the reduction correlated with the locus exhibiting the largest allele frequency fluctuations. These findings highlight fluctuating selection as an important factor shaping genome-wide patterns of genetic variation and effective population size.

The acidic blackwaters of Southeast Asias peat-swamp forests represent some of the most extreme freshwater environments on Earth. Despite their very low pH values, limited nutrients, and hypoxic conditions, these blackwater habitats harbor a remarkable diversity of freshwater fishes, including multiple lineages that have independently adapted to these extreme conditions and, in some cases, exhibiting extreme body miniaturization. These replicate evolutionary lineages therefore provide a powerful comparative framework to investigate adaptation to extreme environments and the genomic basis of miniaturization. Here, we present high-quality, annotated reference genomes for four cypriniform species endemic to these peat-swamp forest ecosystems: Paedocypris sp., Sundadanio atomus, Boraras brigittae, and Rasbora kalochroma. The first two are progenetic miniatures, including Paedocypris, comprising the smallest known fish, while B. brigittae represents a proportioned dwarf and R. kalochroma a non-miniature taxon. Genome sizes ranged from 401-1,290 Mb and heterozygosity from 0.34-1.7%. All genome assemblies achieved pseudo-chromosome-level contiguity, high k-mer completeness (>99%), and high BUSCO completeness (94.5-98.9%). Repeat analyses revealed lineage-specific differences in transposable element landscapes and abundances, while gene annotation identified notable intron length reduction in progenetic miniatures.

Evolution of drug resistance to one drug can alter the minimum inhibitory concentration to another drug. This phenomenon, known as a collateral effect, can manifest as either cross-resistance or collateral sensitivity. Various patterns of collateral effects have been observed experimentally. Repeated adaptation from the same parental strain may result in variable collateral effects; this is non-repeatability. Additionally, adaptation of a pathogen to one drug may produce a specific collateral effect to a second drug, while altering the order of drug exposure can result in a different, or even absent, collateral effect. This phenomenon is termed unidirectionality. The genetic and evolutionary mechanisms underlying these patterns remain incompletely characterised. Here, we propose a frame-work that integrates pharmacodynamics and population genetics and provide minimal examples to explain these patterns and their combinations. Furthermore, we demonstrate that drug concentration and selection regime strongly influence patterns of collateral effects, including repeatability, directionality, and their temporal dynamics.

Nakaseomyces glabratus is a globally distributed opportunistic fungal pathogen. An ongoing discussion in studies of N. glabratus population structure has been whether genetic clusters are best defined using multilocus sequence typing (MLST) or short-read whole-genome sequencing (WGS). To assess the concordance between MLST- and WGS-based phylogenies, we analyzed a dataset of 548 N. glabratus WGS sequences from 12 countries. Clusters identified from WGS largely recapitulated the MLST-defined sequence type (ST) groups: fourteen WGS clusters were composed of a single MLST ST, and the remaining contained STs with very closely related MLST profiles. We thus propose a pragmatic naming convention, consistent with the system used in other microbial species, which specifies WGS cluster labels based on the primary ST. From the large WGS isolate dataset, we determined the prevalence of admixture and genomic variants. Interestingly, seven of the nine singleton isolates were admixed, in addition to 58 isolates from six different clusters. Aneuploidy was detected in 4% of isolates, most commonly in chrE, which contains ERG11, the gene encoding the enzyme targeted by azole antifungals. Aneuploid chromosomes did not exhibit elevated heterozygosity relative to the sequencing error rate, consistent with instability of extra chromosome copies. Copy number variants were found in 3% of the isolates; some of the CNVs co-occurred with aneuploidies, and were primarily identified on chrD, chrE, chrI, and chrM. Our findings demonstrate that deep splits between clusters preserve the utility of MLST ST designations for clade-level designation, yet underscore the utility of WGS for high-resolution genomic analyses. Article SummaryThere is an ongoing debate in studies on Nakaseomyces glabratus about whether traditional MLST analysis is sufficient to determine population structure, or whether the precision of whole genome sequencing (WGS) is necessary. We analyzed WGS data from 548 isolates from around the world. We found a very strong agreement between the two methods. We propose a hybrid naming system, where cluster names are based on the dominant MLST group. We used the WGS data to show that admixed isolates, and those with extra chromosomes or CNVs are rare (<7% of isolates in each class) and are distributed throughout the phylogeny.

Introgressive hybridization between wild and domestic animals is a widespread phenomenon with important implications for genetic diversity, local adaptation, and conservation management. The causes and consequences of this process are poorly understood. In Australia, hybridization between dingoes and domestic dogs presents a dual conservation challenge, threatening the genetic integrity of dingoes while allowing potential adaptive introgression. To investigate the environmental drivers of this process, we analyzed high-density SNP array data in 390 dingoes and 396 domestic dogs. Dingo populations showed regional genetic structure and were clearly differentiated from domestic dogs. Using local ancestry inference and genome-environment association analyses, we found low levels of dog introgression in dingoes from remote areas in Central and Western Australia, and moderate levels in Eastern and Southern populations. Climatic variables (maximum temperature of the warmest month, mean temperature of the driest quarter) and the Human Footprint Index (reflecting density of human populations and environmental modifications) were significant predictors of introgression. We identified four genomic regions with overrepresented dog ancestry, including a large introgressed block on chromosome 27, which contained an olfactory receptor gene showing signatures of positive selection, suggesting adaptive introgression. In addition, a chromosomal inversion previously described in dogs and absent in dingoes was initially identified as an introgressed block. We also detected eight genomic regions nearly free of dog ancestry, suggesting purifying selection against maladaptive variants. Together, these results highlight the complex interplay between introgression, human influence, and local adaptation in dingoes, offering valuable insights for conserving the evolutionary potential of this apex predator in increasingly modified landscapes.

In Toxoplasma gondii, microneme proteins (MICs) are secreted components of the apical complex that play central roles in motility, host cell attachment, and invasion. Because proteins at the host-parasite interface are often predicted to evolve rapidly, MICs have been suggested as candidates for adaptive diversification. We tested this expectation using comparative analyses of three relatively understudied microneme proteins, MIC13, MIC12, and MIC16. Coding sequences were assembled from GenBank and ToxoDB, aligned by translation, and analyzed using maximum-likelihood phylogenetics, codon-based tests of selection, and predicted protein structure. MIC13 was represented by 51 sequences, MIC12 by 30, and MIC16 by 34, spanning multiple T. gondii haplogroups and including Hammondia hammondi and Neospora caninum as outgroups. All three genes were highly conserved among T. gondii strains, but their phylogenetic trees were topologically incongruent, indicating that individual MICs do not recover a single shared strain history. Contrary to expectation, no positively selected codons were detected in any gene. Instead, purifying selection was detected at one site in MIC13 and 15 sites in MIC12, while no significant codon-specific selection was detected in MIC16. Several constrained MIC12 sites overlapped annotated EGF and calcium-binding EGF-like domains, consistent with structural conservation of extracellular adhesion modules. AlphaFold prediction of MIC13 supported two sialic acid-binding micronemal adhesive repeat regions, but the single constrained MIC13 site did not overlap these motifs. Together, these results indicate that MIC13, MIC12, and MIC16 are shaped more by sequence conservation and heterogeneous gene histories than by strong recurrent positive selection. These findings refine expectations for microneme evolution in T. gondii and highlight conserved domains that may be important for parasite invasion and future functional study.

Plasmids are extrachromosomal mobile genetic elements that can facilitate rapid bacterial adaptation by transferring genes between individuals. While plasmids are known to exist in diverse habitats and encode a range of traits, most of our knowledge about plasmids comes from clinically-associated antimicrobial resistance (AMR) plasmids that have already been recruited as vectors of drug resistance and have likely been shaped by strong selection for plasmid-encoded resistance. Here, we investigated 26 plasmids from the pQBR collection -- a set of large, co-existing mercury resistance environmental plasmids isolated in Pseudomonas spp. from a field in Oxfordshire in the 1990s -- and explored the ability of pQBR plasmids to mobilise novel chromosomally-encoded traits. New whole genome sequences for 25 plasmids confirmed that these soil-isolated plasmids are generally very large (140-588 kb), constitute at least five distinct genetic groups, and have relatives in various other Pseudomonas species and habitats. Despite significant nucleotide-level divergence, Groups I (pQBR103-like, [~]406 kb) and IV (pQBR57-like, [~]328 kb) showed remarkable ancient similarities in synteny and gene content both with one other, and with the PInc-2 family of plasmids known to mobilise clinically significant drug resistance in Pseudomonas aeruginosa. None of the pQBR plasmids sequenced to date harboured known AMR determinants, but putative phage defence systems and metal resistances were evident. Transposable elements, including the Tn5042 mercury resistance transposon, were responsible for significant structural variation within plasmid groups, consistent with a predominant role of transposons in rapidly remodelling plasmids. To experimentally test the ability of pQBR plasmids to spread new traits, we developed a novel transposon mobilisation assay which showed that certain Group IV pQBR plasmids were especially effective at acquiring the chromosomally-encoded transposon Tn6291, and that this mobilisation was likely due to specific plasmid factors rather than generic conjugation rate. Our work presents a tractable set of sequenced plasmids suitable for exploring the evolution and dynamics of gene acquisition by pre-AMR plasmids, and provides a key case study highlighting the pervasive interplay between plasmids and transposable elements that can drive microbial genome evolution. Repositories: github.com/jpjh/PQBR_PLASMIDS Impact statementPlasmids can drive microbial evolution by acting as vectors for horizontal gene transfer. Because of their central role in disseminating antimicrobial resistance (AMR), plasmids are mainly explored as vehicles for AMR traits, meaning that our knowledge of the diversity and evolutionary dynamics of non-AMR plasmids is more limited. Here, we explore sequences from a set of mercury resistance plasmids isolated in Pseudomonas spp. from pristine agricultural land that lack AMR determinants. By providing new whole genome sequencing analyses we expand the set of sequenced pQBR plasmids to 26, finding globally dispersed relatives from clinical, environmental, and industrial settings, and identifying an ancient plasmid backbone shared amongst divergent modern environmental and clinical AMR plasmids. We experimentally verify the role of pQBR plasmids in readily mobilising chromosomal traits using a novel transposon mobilisation assay, which suggests that specific plasmid-transposon interactions may drive trait spread. Overall, our work expands our understanding of the role of environmental plasmids in mobilising and disseminating adaptive traits.

Transitions from sexual to asexual reproduction are well-documented across different taxa. However, despite extensive efforts, the regulatory changes underlying the emergence of asexuality remain largely undiscovered in the majority of species studied. Artemia brine shrimp have multiple closely related sexual and obligate parthenogenetic lineages, making them a promising model for addressing this question. While earlier work suggested that asexuals use a modified meiosis, and inferred a likely role for the Z-chromosome in its transmission, no master regulator or genetic changes have been put forward as the root causes for the shift. Here, we generate single-nucleus RNAseq data of the female reproductive system of individuals from the Aibi lake population of Artemia parthenogenetica and its closely related obligate sexual species Artemia sp. Kazakhstan. We identify the germline cell clusters in the female reproductive system and perform differential expression analysis to infer substantial transcriptional differences at genes putatively involved in cell cycle and oocyte development between the meiotic cells of the two species. Additionally, we use whole-genome sequencing of 32 individuals from two backcrossing experiments to narrow down the genomic regions associated with the transmission of asexuality to an 8 megabase region of the Z chromosome. Within the identified regions, two adjacent genes with known functions in oogenesis, ITPR and USP8, show differential expression and genetic differentiation between sexuals and asexuals, making them promising candidate drivers of asexuality in this species. Significance statementWhile most animals reproduce sexually, many do not, and why and how these shifts occur remains an open question. This paper presents a systematic investigation of the molecular changes that underlie the transition from sexual to asexual reproduction in brine shrimp. We combine multiple computational and experimental approaches to look for differences between close sexual and asexual lineages. We find that a subset of meiotic germ cells is regulated differently in the two, and that two important oogenesis genes are the likely drivers of asexuality. This work is unique in providing an in-depth characterization of the combined genetic and regulatory changes underlying this key transition in reproductive modes.

It has been known for decades that bacteriophages encode tRNA genes, but their function and the factors contributing to their acquisition and retention are unclear. Although tRNAs are found in a variety of phages infecting a variety of bacteria, many large-scale computational studies investigating tRNA acquisition and retention in phages are specific to Mycobacterium phages; however, these findings may not be representative of other phages or bacteria. This work uses a broader sampling of phages and hosts to investigate the relationships between codon usage bias, infection cycle, and tRNA gene numbers in phage genomes. We analyzed 154 phages infecting 7 host genera, including Gram-negative (Escherichia, Shigella, Salmonella) and Gram-positive (Bacillus, Lactobacillus, Staphylococcus, Mycobacterium) bacteria. Phages included temperate and virulent representatives, plus a range of tRNA numbers and morphologies. All phages and hosts were analyzed using four metrics: GC content, Effective Number of Codons, Relative Synonymous Codon Usage, and tRNA Adaptation Index. On a global scale, virulent phages with many tRNA genes show greater differences in codon usage and codon adaptation compared to their respective hosts. Gram-negative bacteria and their phages generally exhibit greater differences in codon usage compared to Gram-positive bacteria and their phages. Phages infecting Gram-negative hosts also tend to encode more tRNA genes. In nearly all genus-level comparisons, Mycobacterium phages were different from any other host and from global patterns. This suggests previous computational studies performed in Mycobacterium phages are likely not applicable on a global scale or to phages infecting other host genera. AUTHOR SUMMARYBacteriophages, or phages, are viruses infecting bacteria. They are abundant in all environments, yet how they interact with their bacterial hosts is still not well-understood. Like other viruses, phages must rely on the host translational components to replicate and form new phage particles; and similarly to other parasites, phages have genomes that differ significantly from their hosts in terms of composition. In this work, we explore the relationship between phage lifestyle, number of tRNA genes encoded, and genome differences from the host using a variety of phages and their associated hosts. Phages can be either virulent (do not integrate into the host genome) or temperate (capable of integrating into the host genome), with differences from the host genome more pronounced in virulent phages. There are many phages that also carry tRNA genes, and having higher numbers of tRNAs is associated with larger differences from the host genome. The findings here indicate that virulent phages carrying large numbers of tRNAs diverge the most from host genome composition.

The loss of an organelle represents an extreme case of reductive evolution. Only a few plastid loss events have been documented and they seem to be particularly rare in free-living ancestrally photosynthetic taxa. We isolated a new phagotrophic protist, Ciliophrys sp. Baltic, representing a poorly studied non-photosynthetic lineage in the stramenopile algal class Dictyochophyceae. Genomic and transcriptomic data obtained from this organism enabled us to demonstrate that it lacks a plastid genome and nuclear genes critical for plastid biogenesis. The plastid organelle is thus most likely secondarily missing from Ciliophrys sp. Baltic, although a few enzymes originally located in the plastid persist with a changed localization. Unexpectedly, Ciliophrys sp. Baltic harbors a bacterial endosymbiont, Candidatus Penulousia baltica sp. nov. (Rickettsiales). Analysis of its complete genome sequence uncovered unique aspects of its metabolism not previously reported in other Rickettsiales. The bacterium exploits its host metabolically and presumably manipulates it with an arsenal of putative effector proteins. However, we also found that it may help offset some of the biosynthetic deficiencies of its host stemming from the plastid loss by providing haem and a lysine precursor. Altogether, our study establishes a new single-cell model system for studying plastid loss and endosymbiont acquisition.

Organellar genomes are both a resource for reconstructing organismal phylogenies and interesting subjects for evolutionary studies. Herein, we focused on the organellar genomes of eustigmatophytes (eustigs), a class of the algal phylum Ochrophyta with a growing biotechnological potential, and massively expanded the existing limited sample by 51 new organellar genomes. Analyses of this large dataset provided a robustly resolved eustig phylogeny and important insights into the evolution of unique features of eustig organellar genomes. Eustig plastomes are rather stable in terms of the gene content and order, with only minor differences stemming from differential gene loss and rare lineage-specific gain. In contrast, eustig mitogenomes vary broadly in their gene order, the content of "accessory" genes, and substitution rates. Notably, the new data illuminated the origin of some of the organellar genes previously deemed eustig-specific. Thus, the plastid gene ycf95 most likely is an extremely divergent version of ycf35, and the mitochondrial genes orfX and orfY evolved, respectively, by rps4 duplication and extreme divergence of rps1 ortholog. Most interestingly, we identified five previously unrecognized orthogroups of mysterious mitochondrial orfs that are patchily distributed across eustigs yet likely evolved in the ancestor of this class. These orfs have no discernible homologs outside eustigmatophytes but are predicted to encode multipass membrane proteins with a soluble C-terminal domain. Finally, our results revise some of the previous conclusions regarding the mitochondrial translation in eustigs and suggest the recruitment of a group of unusual tRNAs for a translation-independent function in the genus Vischeria.

High and low tides occur twice a day (every [~]12.4 hours), with the largest tidal ranges during spring tides around new and full moons (every [~]14.765 days). While these lunar cycles are known to influence many animal phenotypes, particularly the reproduction of coastal animals, the genetic basis of lunar-related rhythms remains unclear. Since phenotypic variation is a valuable resource for elucidating such mechanisms, we examined geographic variation in the lunar-regulated mass spawning of the grass puffer (Takifugu alboplumbeus) along the Japanese coast. We found that western populations spawn during the first half of the spring tides, whereas eastern populations spawn during the second half. Furthermore, although spawning typically occurs a few hours before high tide, this timing is restricted to a specific time window that is earlier in the western populations than in the eastern ones. Behavioral analysis of larvae also revealed a shorter free-running circadian period ({tau}) in the western population than in the eastern ones. As differences in {tau} affect individual variation in the timing of physiological functions and behaviors, we hypothesized that differences in {tau} could account for the different time windows and consequently the observed difference in spawning days. Population genomics analysis identified proline-rich transmembrane protein 1-like (prrt1l) as a candidate gene. Expression of prrt1l was observed in the circadian pacemaker suprachiasmatic nucleus, and triple CRISPR F0 knockout of prrt1l shortened the free-running period in larvae. These findings suggest a potential mechanism underlying the geographic variation in lunar-synchronized spawning behavior. HighlightsO_LIThe geographic variation exists in the lunar-regulated spawning of the grass puffer, with differences in spawning dates and times between western and eastern Japan. C_LIO_LIThe free-running period of western populations is shorter than that of eastern populations, which is consistent with their earlier spawning timing. C_LIO_LIPopulation genomics analysis identified prrt1l as a candidate gene harboring population-specific missense mutations, the knockout of which shortens the free-running period. C_LI

Parthenogenesis is relatively rare and often regarded as an evolutionary dead end. Despite this, certain parthenogenetic animal species have endured for millions of years, but it is unclear what enables the persistence of some parthenogenetic lineages. Transitions from sexual to parthenogenetic reproduction can occur through different evolutionary processes that give rise to diverse cytological reproductive mechanisms. These mechanisms are likely to influence genetic diversity, especially in the early stages after the transition to parthenogenesis and may thus affect lineage persistence. To understand such evolutionary transitions, we used experimental crosses to investigate the mechanism of parthenogenesis and the immediate genetic consequences of switching from sexual to parthenogenetic reproduction in the facultatively parthenogenetic phasmid Megacrania batesii. We obtained DNA sequence data from multiple lineages propagated over three generations via sex, parthenogenesis, or transitions between reproductive modes. We quantified heterozygosity and within-family genetic variation and compared the genetic patterns with predictions for known mechanisms of parthenogenesis. We found that a single generation of parthenogenesis typically resulted in (near-)complete loss of heterozygosity and an absence of within-family genetic variation, consistent with automixis with gamete duplication or terminal fusion and little/no recombination. However, we also found evidence of variation in the mechanism of parthenogenesis among lineages and even within the same individual, associated with drastic differences in the amount of heterozygosity and within-family genetic variation maintained across generations. Our findings show that considerable variation in parthenogenetic mechanisms can exist within populations and suggest that such variation could influence the persistence and evolution of parthenogenetic lineages.